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ABSTRACT 

We have explored the structure of hot magnetized accretion flow with thermal conduction. The 
importance of thermal conduction in hot accretion flows has been confirmed by observations of the 
hot gas surrounding Sgr A* and a few other nearby galactic nuclei. For a steady state structure of 
such accretion flows a set of self similar solutions are presented. In this paper, we have actually tried 
to re-check the solution presented by Abbassi et al. (2008) using a physical constrain. In this study we 
find that Eq 29 places a new constrain that limits answers presented by Abbassi et al. 2008. In that 
paper the parameter space in which it is established in the new constrain was plotted. However, the 
new requirement makes up only a small parameter space with physically acceptable solutions. And 
now in this manuscript we have followed the idea with more effort, and tried to find out how thermal 
conduction influences the structur of the disks in a physical parameter space. We have found out 
that the existence of thermal conduction will lead to reduction of accretion and radial and azimuthal 
velocities as well as the vertical thickness of the disk, which is slightly reduced. Moreover, the surface 
density of the disk will increase when the thermal conduction becomes important in the hot magnetized 
flow. 

Subject headings: accretion flow, magnetic field, thermal conduction 



1. INTRODUCTION 

It is widely believed that many astrophysical objects 
are powered by mass accretion on to black holes. Black 
hole accretion disk models have been extensively stud- 
ied during the past three decades (see Kato et al. 2008 
for a review). The standard geometrically thin, optically 
thick accretion disc model can be successfully explained 
by most of the observational features in active galactic 
nuclei (AGN) and X-ray binaries (Shakura & Sunyaev 
1973). Besides the traditional standard disk model by 
Shakura and Sunayev (1973), there are a number of new 
disk classifications including advection-dominated accre- 
tion flows (ADAF) or radiative inefficient accretion flows 
(RIAF) for a very small mass-accretion disks (Narayan & 
Yi 1994) and supercritical accretion disk or so-called slim 
disk, for a very large mass- accretion rate (Abramowicz et 
al 1998). ADAFs are optically thin and hot (compared 
to the virial temperature of the gas in the flows), which 
radiate mostly in X-ray waveband (see Narayan & Mc- 
Clintock 2008). In the ADAF model, only a small frac- 
tion of the gravitational energy released in the accretion 
flow is radiated away due to inefficient cooling and most 
of it energy is stored in the accretion flow and advected 
to the black hole. This model can successfully explain 
the main observational features of black hole X-ray bina- 
ries and low- luminosity AGN (LLAGN) (e.g., Narayan 
& Yi 1994, 1995a; Quataert et al. 1999; Yuan et al. 
2003; Ho 2008). The diversity of models indicates that 
modeling the hot accretion flows around a black hole is a 
challenging and controversial problem. One of the largely 
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neglected physical phenomena in the modeling ADAFs is 
the thermal conduction. Recent observations of hot ac- 
cretion flows around active galactic nuclei indicate that 
it should be based on collision-less regime (Tanaka & 
Mcnou 2006) . It is also suggested that accretion in these 
systems proceeds under the weakly collisional condition. 
Furthermore, it is proposed that thermal conduction can 
be a possible mechanism through which sufficient extra 
heating is provided in hot advection dominated accretion 
flows. So, thermal conduction seems to have an impor- 
tant role in energy transport in the accreting materials 
around a hot accretion disk where they are nearly com- 
pletely ionized. Since advection-dominated disks have 
high temperatures, the internal energy per particle is 
high. This is one of the reasons why advective cooling 
overcomes radiative cooling. For the same reason, tur- 
bulent heat transport by conduction in the radial direc- 
tion is non-negligible in the heat balance equation. So it 
should be important to consider the role of thermal con- 
duction in an ADAF solution. Shadmehri (2008), Faghei 
(2011), Abbassi et al. (2008, 2010) Tanaka & Menou 
(2006) have studied the effect of hot accretion flow with 
thermal conduction in a semi- analytical method; physics 
of such systems have been studied in simulation models 
(e. g. Sharma et al. 2008; Wu et al. (2010). Shadmehri 
(2008) has shown that thermal conduction opposes the 
rotational velocity, but increases the temperature. Ab- 
bassi et al. (2008) have shown that for this problem there 
are two types of solutions; high and low accretion rates. 
They plotted the radial velocity for both solutions which 
revealed that it is influenced by thermal conduction. In 
this study, by introducing a new physical condition, it has 
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been shown that high accretion rate solutions of Abbassi 
et al. (2008) are not exactly correct. Also some of the low 
accretion rate solutions presented there are not of physi- 
cal value. With this extra condition it is possible to find 
the physically meaningful parametric space using which 
we next plotted the dynamical quantities to investigate 
how thermal conduction affects them. Finally, tangled 
magnetic fields in accretion flows are likely to reduce the 
effective mean free path of particles. The magnitude of 
this reduction, which depends on field geometry, is still 
unknown. In this paper, we will investigate the effect of 
thermal conduction on the physical structure of ADAF- 
like accretion flow around a black hole in the presence of 
a toroidal magnetic field. 

2. THE BASIC EQUATIONS 

We are interested in analyzing the structure of a mag- 
netized ADAF where thermal conduction plays an im- 
portant role in energy transportation. We only consider 
toroidal component of the magnetic field and suppose 
that the gaseous disk is rotating around a compact object 
of mass M*. Thus, for a steady axi-symmetric accretion 
flow, i.e. d/dt = d/d(f> = 0, we can write the standard 
equations in the cylindrical coordinates (r, <j>, z). In ad- 
dition, we vertically integrate the flow equations and so, 
all the physical variables only become functions of radial 
distance r. Moreover, relativistic effects are neglected 
and Newtonian gravity in radial direction is taken into 
account. The disk is supposed to be turbulent and pos- 
sesses an effective turbulent viscosity v. We adopt a- 
prescription for viscosity of a rotating gaseous disk. As 
for conservation of energy, it is assumed the energy gen- 
erated due to viscosity dissipation is balanced by the ra- 
diation, thermal conduction and advection cooling (e.g., 
Narayan & Yi 1994, Abbassi et al. 2008, 2010). The 
equation of continuity reads 
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where V r is the accretion velocity (V r < 0) and S = 2pH 
is the surface density at a cylindrical radius r. Where p 
the mass loss rate per unit volume, H is the disk's half- 
thickness. The equation of motion in the radial direction 
is: 
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where v v , c s and ca are the rotational velocity of the disk, 
sound velocity and Alfven velocities of the gas, respec- 
tively. Sound speed and the Alfven velocity arc defined 
as c 2 = p gas /p and c\ = B^/Airp = 2p mag /p where p gas 
and p mag are the gas and the magnetic pressures, respec- 
tively. 

The vertically integrated angular momentum equation 
becomes: 

where a is the viscose parameter, 17 (= ^) the angular 
velocity and fl k the Keplerian angular velocity, respec- 
tively. 

By integrating over z of the hydrostatic balance one 
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where j3 = ^"" g = 5 (7-) 2 which indicates the impor- 
tant of magnetic field pressure compared to gas pressure 
(Akizuki & Fukue 2006, Abbassi et al. 2008, 2010). Two 
cases can be considered; Case 1: where the pressure is as- 
sumed to be the gas pressure (thermal pressure) (which 
is the choice of this paper) and Case 2: when the pressure 
is assumed to be the total of magnetic plus gas pressure. 
So in case 2 we replace a with a(l + /3) in all of the equa- 
tions. Next, the dynamical properties of the disk for dif- 
ferent values of /3 is demonstrated. The energy equation 
considering cooling and heating processes in an ADAF 
can be found as follows. We assume the energy gener- 
ated due to viscous dissipation and the heat conducted 
into the volume are balanced by the advection cooling 
and energy transport by thermal conduction. Thus, 

7-1 dr 2IiVr ° s dr~ n k r( dr> r dr (r * s) 

(5) 

On the right-hand side of the energy equation we have: 

Q+ Q— Qcond Qf "f" Qcond 

where Q con d = — V • F s represents energy transfer due 
to the thermal conduction and F s = 5& s pcl is the sat- 
urated conduction flux ( Cowie & Mackeel977) in the 
direction of temperature gradient. 4> s , the dimensionless 
coefficient, is less than unity. Finally since we consider 
the toroidal component for the global magnetic field, the 
induction equation can be written as: 
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where B v is the field scaping/creating rate due to mag- 
netic instability or dynamo effect. Now we have a set 
of MHD equations that describe the structure of mag- 
netized ADAFs. The solutions to these equations are 
strongly correlated to viscosity, magnetic field strength 
/?, the degree of advection / and the thermal conduction 
parameter 4> s for the disks. We seek a self-similar solu- 
tion for the above equations. In the next section we will 
present self-similar solutions to these equations. 

3. SELF-SIMILAR SOLUTIONS 

In order to have a better understanding of the physical 
processes taking place in our disks, we seek self-similar 
solutions of the above equations. The self-similar method 
has a wide range of applications for the full set of MHD 
equations although it is unable to describe the global 
behavior of accretion flows since no boundary conditions 
have been taken into account. However, as long as we 
are not interested in the behavior of the flow near the 
boundaries, these solutions are still valid. In the self- 
similar model the velocities are assumed to be expressed 
as follows: 

V r (r) =- Cl aV k (r) (7) 
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and constant ci,c 2 and c 3 are determined later from the 
main MHD equations. We obtain the disk half-thickness 
H as: 

— = a/c 3 (1 + /3) = tana (12) 

Hence, a hot accretion also has a conical surface, whose 
opening (half-thickness) angle is a. Assuming the surface 
density E to be in the form of: 



we obtained, e.g. 
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By substituting the above self-similar solutions in to 
the the equations of the disks, we obtain the following 
system of dimensionless equations, to be solved for C\,C2 
and c 3 : In case 1: 
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It is evident from the above equations that for s = 
— 1/2 there is no mass loss while it is present for the 
case where s > —1/2. On the other hands, the escape 
and creation of magnetic fields balance one another for 
s = 3. Although outflow is one of the most important 
processes in accretion theory,(see Narayan & Yi 1995; 
Blandford & Begelman 1999; Stone & Pringle & Bcgcl- 
manl999 and also some recent work like Xie & Yua 2008; 
Ohsuga & Mineshige 2011), but in our model we choose 
a self-similar solution in which p — and B$ oc r~ n / 4 
(s = —1/2) thus ignoring the effect of wind and outflow 
on the structure of the disk. After some algebraic ma- 
nipulations, we obtain a forth order algebraic equation 
for c\\ 

D 2 c\ + 2DBc\ + {B 2 -2D)c\-{2B + A 2 ) Cl + l = (21) 

where the coefficients depend on the input parameter 
as: 



D = -a 2 
2 



(22) 



9/7-1 2' 



3(s + l) J 



9fa y 2 ;L 3(s + l) J 



(23) 
(24) 



This algebraic equation shows that the variable C\ 
which determines the behavior of radial velocity only de- 
pends on the a, $ s , (3 and /. Using c\ from this equation, 
the other variables (i.e. c 2 and c 3 ) can be determined 
easily: 
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As was stated above we are solving the problem in the 
case where there is no wind, s = —1/2. When s = — \ 
the expressions for c\ and c 2 are found to be: 
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Abbassi ct al. (2008) have found two types of solutions 
which represent high and low accretion rates. These two 
rates act differently towards the saturated thermal con- 
duction parameter, <j> s . As it was stated in the introduc- 
tion the main aim of this paper is to place an additional 
physical condition. This new requirement limits the so- 
lutions of c 2 > 0. From the above equation, c 2 > and 
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|, we have: 
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Therefore the (cj) s , a) parametric space could be found as 
follows: 
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A numerical solution to this equation is presented in Fig. 
1 for variety of values of the viscosity parameter a and 
the thermal conduction parameter <f> 8 . These solutions 
may be considered either as flows with a fixed value of 
viscosity parameter and a sequence of increasing thermal 
conduction parameter or a sequence of different values 
of viscosity parameter a in a fixed thermal conduction 
regime. These plots show that the new physical con- 
dition places a constraint on the physically valid para- 
metric space. As a result we used this extra condition 
to check the validity of the solutions presented by Ab- 
bassi et al. (2008). Tanaka & Meneu (2006), Abbassi 
ct al. (2008) have shown that for a very small <j) s their 
solutions coincide the original ID ADAF solutions; but 
by adding the saturated conduction parameter, <fr s , the 
effect of thermal conduction can be better seen while ap- 
proaching 0.001 — 0.01. On the other hand, the widely 
accepted values of a are between 0.01 — 0.1. So we have 
to plot our solutions in this range. The dark areas in Figl 
are where we are looking for; regions of the parametric 
space physically acceptable. Another aim of this study 
was tore-check the solution presented by Abbassi et al. 
2008 and and further to study the effect of heat conduc- 
tion and magnetic field on the structure of ADAFs. Fig. 
1 shows that for larger values a we should choose a large 
<f> a in order to be sure to have a physically valid solution. 
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Fig. 1. — The parameter space for thermal conduction parameter 
4> s as a function of viscous parameter a. The dark section of the 
plot shows the physically accepted area of the parametric space 
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Fig. 2. — Numerical coefficient ds as a function of advection 
parameter / for several values of a viscose parameter. All of this 
figures was set up for a = 0.07 ' ,<f> s = 0.003 

The calculations are carried out in a range of viscosities, 
thermal conduction and magnetic field parameters. The 
calculated global structures of some ADAFs with differ- 
ent parameters are plotted in Figures 2-8. 

Figure 2 shows how the coefficient C,s depend on the 
magnetic field parameter, (3 for different values of advec- 
tion parameter /. Radial velocity is determined by Cf 
which is shown in the upper panel. In ADAFs the ra- 
dial velocity is generally less than free fall velocity on a 
point mass, but it becomes larger if the advection param- 
eter / is increased and such behavior is consistent with 
the previous analytical solutions (e.g., Akizuki & Fukue 
2006; Abbassi et al. 2008, 2010). For a given /, both 
the radial and rotational velocities increase as the disk 
becomes more magnetized (i.e., larger /3). A magnetized 
disk must rotate faster than the case where no magnetic 
field is present because of the effect of magnetic tension. 
It is clear that increasing the magnetic field strength in- 
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Fig. 3. — Numerical coefficient ds as a function of advection 
parameter / for several values of 4> s thermal conduction parame- 
ter. All of this figures was set up for a = 0.07,/3 = 0.1 
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Fig. 4. — Numerical coefficient ds as a function of thermal con- 
duction parameter ij>s for several values of /3 magnetic field param- 
eter All of this figures was set up for a = 0.07,/ = 0.7 

creases radial velocity (accretion rate) in the accretion 
flows. Radial velocity of the disk is more sensitive to 
the variations of / when the disk is more magnetized. 
The rotational velocity behaves the same for different 
values of / but it slightly shifts up when parameter /? in- 
creases. The standard ADAFs becomes thicker when the 
disks is advective / ~ 1 but by adding the magnetic field 
strength the vertical thickness is increased significantly 
for a given /. In Fig. 3 we investigate the role of satu- 
rated thermal conduction in radial, toroidal and the ver- 
tical thickness of accretion flows. Increasing the thermal 
conduction coefficient <j) s will decrease the radial velocity. 
It will also have a large effect on the rotational velocity 
of accretion flows (middle panel). Furthermore it will de- 
crease the sound speed and therefore vertical thickness of 
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Fig. 5. — Numerical coefficient ds as a function of thermal con- 
duction parameter <j> a for several values of / advection parame- 
ter. All of this figures was set up for a = 0.07,/9 = 0.1 

the disks too (lower panel) . In Figs 4 and 5 we have plot- 
ted the coefficient Ci with thermal conduction parameter 
<p s for different values of /3 and f. Tanaka &Menou (2006) 
have shown that for a very small <p s their solutions co- 
incide the original ADAF solution (Narayan & Yi 1994), 
but by adding the saturated conduction parameter, <j> 3 , 
the effect of thermal conduction can be better seen while 
approaching 0.001 — 0.01. On the other hand, the widely 
accepted values of a are between 0.01 — 0.1. So we have 
to plot our solutions in this range. We can see that for a 
given set of input parameters, the solution reaches a non- 
rotating limit at a specific value of (j> s . We cannot extend 
the solution for values of this limit for (f> s because equa- 
tion (25) gives a negative value for Cf , which is clearly 
not valid. The values of this limitation on (j) s increases by 
adding magnetic field parameter, (3, fig 4 (middle panel), 
advection parameter, /, fig 5 (middle panel) smoothly. 
In fact the higher values of f3, / will allow the disk to 
have a physical rotating solution for larger <fi s . As it is 
seen in the lower panels of figures 4 and 5 by increasing 
/3, / the vertical thickness of the disks in a given (£> s will 
increase. But the vertical thickness in a fixed /? or /, will 
decrease globally as <fr s increased. It would be expected 
that in high <f> s , disks becomes thin as it was predicted 
by Abbassi et al (2010) and Shadmehri (2008). 

As the level of thermal conduction is increased, more 
heat flows outward from the hot inner regions, resulting 
in a local increase of the gas temperature relative to the 
original ADAF solutions. Simultaneously, the gas ad- 
justs its angular velocity and increases its inflow speed 
to conserve its momentum balance. We have found that 
the level of advection, /, will significantly modify the 
structure of the disk. The breakdown of the solutions, 
when Cf — > 0, occurs at lower values of (f) s . Clearly 
thermal conduction can significantly affect the structure 
and the properties of the hot accretion flows. In order 
to understand the role of thermal conduction better it 
is useful to study it in a more realistic two-dimensional 
case, without height integrations. But at this stage we 
are going to find out how thermal conduction affects the 
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FlG. 6. — Surface density T,(kg/m 2 ) as a function of the radius 
(r)of disk for several values of <p s , f, /3, a parameters 





a=0.09 
i a=0.08 


f=0.7 , p=0.1 , ♦ =0.003 












stV 




r 


10 




f=0-7 


a=0.07 , (3=0.1 , + s =0.003 












o' 




r 


10 




S =0.004 <i s =0.003 


a=0.07 ,p=0.1 ,f=0.7 
















r 


10 




, P=0.5 


a=0.07,f=0.7,+ =0.003 







Fig. 7. — Radial velocity v r (m/s) as a function of the radius (r)of 
disk(in unit of r g )for several values of <p B , f, /3, a parameters 

structure of the disks with a simple analysis. It would be 
interesting to show how the surface density, radial veloc- 
ity and accretion rate which are all observables, change 
with the radius of the disk when the thermal conduction 
plays an important role. To show how surface density is 
affected by our input parameters we use: 

M = 27rr£(-7j r ) (31) 

where accretion rate is constant from the continuity 
equation; Eq(l). Assuming M w 10 _8 M Q yr _1 and 
M = 10M Q and by using self-similar solution for v r we 
have: 

S = r-i (32) 

2TrC\aVGM 

As it is seen in figures 6 and 7, we have plotted the sur- 
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face density, E, and radial velocity as a function of r. 
As expected, the absolute values of surface density and 
radial velocity will decrease as one moves out farther in 
the disk. In figure 6 it is evident that surface density 
is increased when thermal conduction plays an impor- 
tant role. But by adding the advection parameter £ will 
decrease significantly. By adding viscous and advection 
parameters the absolute values of radial velocity will in- 
creased in a fixed r, as expected. Thermal conduction 
has the same effect and causes high radial velocity when 
it has an important role. Compared to other physical 
parameters, the magnetic field, /3, does not have a signif- 
icant effect on surface density and velocity gradients. 



4. SUMMERY AND CONCLUSION 

In this paper, we studied an accretion disk in the advec- 
tion dominated regime by considering a purely toroidal 
magnetic field and in presence of thermal conduction. 
Some approximations were made in order to simplify 
the main equations. We assumed an axially symmetric, 
static disk with the a-prescription of viscosity. We also 
ignored relativistic effects and self-gravity of the disk. 
Considering the weakly collisional nature of a hot accre- 
tion flow (Tanaka & Menou 2006; Abbassi et al. 2008), a 
saturated form of thermal conduction was adopted as a 
possible mechanism for energy transportation. We have 
accounted for this possibility by allowing the saturated 
thermal conduction constant, <p s , to vary in our solutions. 
A set of similarity solutions was presented for such a con- 
figuration. Also following those presented by Abbassi 
et al. (2008) we re-examined the validity of their solu- 
tions by putting an additional physical condition. They 
had found two types of solutions: low and high accretion 
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rate solutions. Their solutions had different behaviors 
towards thermal conduction parameters. We found out 
that the high accretion rate solution presented by Ab- 
bassi et al. (2008) is not physically valid. Also our new 
imposed condition puts some physical constrain on low 
accretion rate solutions. As a result in order to have a 
physically valid solution only a small part of parameter 
space should be considered. After we assuring the va- 
lidity of our solutions, by using proper values of input 
parameter, a, (3, (j) s , /, we investigated the influence of 
thermal conduction on the structure of accretion flows. 
It is shown that the surface density is increased when the 
thermal conduction plays an important role. Also ther- 
mal conduction has the same effect on radial velocity. 
However, there exist some limitations to our solutions. 
One is that the self-similar hot accretion flow with con- 
duction is a single-temperature structure. Thus if one 
uses a two-temperatures structure for the ions and elec- 
trons in the disks, it is expected that the ions and elec- 
tron temperatures decouple in the inner regions, which 
will modify the role of conduction. The other limitation 
of our solution is the anisotropic character of conduction 
in the presence of magnetic field. Balbus (2001) has ar- 
gued that the structure of the hot flows could be affected 
by the anisotropic character of thermal conduction in 
the presence of magnetic field. Although our preliminary 
self-similar solutions are too simplified, they clearly im- 
prove our understanding of the physics of ADAFs around 
a black hole. To have a realistic picture of an accre- 
tion flow a global solution is needed rather than the self- 
similar one. In our future studies we intend to investigate 
the effect of thermal conduction on the observational ap- 
pearance and properties of a hot magnetized flow. 
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